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CL|- ABSTRACT 

< : 

Aims. We have selected the Galactic Hn region M43, a close-by apparently spherical nebula ionized by a single star (HD37061, 
C . B0.5 V) to investigate several topics of recent interest in the field of H n regions and massive stars. In a series of two papers, we 

perform a combined, comprehensive study of the nebula and its ionizing star by using as many observational constraints as possible. 

Methods. We collected for this study a set of high-quality observations, including the optical spectrum of HD 37061, along with 

nebular optical imaging and long-slit spatially resolved spectroscopy. The first part of our study comprises a quantitative spectroscopic 
r ^ , analysis of the ionizing star, and the empirical analysis of the nebular images and spectroscopy. All the information obtained here will 

■ be used to construct a customised photoionization model of the nebula in Paper II. 

' Results. We determine the stellar parameters of HD 37061 and the total number of ionizing photons emitted by the star. We find 

O l| observational evidence of scattered light from the Huygens region (the brightest part of the Orion nebula) in the M 43 region. We 

show the importance of an adequate correction of this scattered light in both the imagery and spectroscopic observations of M 43 
in accurately determining the total nebular Ha- luminosity, the nebular physical conditions, and chemical abundances. We perform a 

' detailed nebular empirical analysis of nine apertures extracted from a long-slit located to the west of HD 37061 (east-west direction), 

obtaining the spatial distribution of the physical conditions and ionic abundances. For three of the analyzed elements (O, S, and N), 
we determine total abundances directly from observable ions (no ionization correction factors are needed). The comparison of these 
abundances with those derived from the spectrum of the Orion nebula indicates the importance of the atomic data and, specially in 

' the case of M 42, the considered ionization correction factors. 

Key words. ISM: HII regions - ISM: individual: M43 - ISM: abundances - Stars: early-type - Stars: fundamental parameters - 
Stars: atmospheres - Stars: individual: HD 37061 - 

: 

m . 

1. Introduction of the primary, they can be neglected in terms of ionization of 

' the nebula. 

U M43 (NGC1982) is an apparently spherical Hn region sev- M4 3 belongs to t he same molecular complex as the Orion 

" . eral arcmtn to the northeast of the well-known Orion nebula. nebula riGoudislll982l) . In particular, M 43 is located at the north- 

HD 37061 (NU On, Par 2074, Brun 747), an early B-type star east border of the extended Orion nebula (EON), a large ellip- 

wtth broad lines located at the center of the nebula, is the main tical stmcture surround ing the Huygens region (corresponding 

ionizing source of M43. The spectral classifications found in the to the centralj brighter par t of the Orion nebula). Although the 

literature for HD 37061 range between B0.5 V and Bl V, prob- p rec ise structure of th e EON is still un known, recent studies 

,3 i ably because of the difficulty in detecting the faint He n4541 to'Dell & GossL 1207)1 lO'Dell & Harrisl |S) have associated 

line when using low resolution spectra and photographic plates, the dillusd]] and extended emission arising from this region with 

Up to three componentshave been identifiedwithin this stel- scattere d light, produced by dust reflecting the stellar contin- 

lar system ( |Abt et alj [199J | Preibisch et aU | 1999| ). The primary uum emission from the Trapezium cluster stars and the nebular 

star is known to be a spectroscopic binary (SBD with an esti- line emission from the Huygens region. A dust lane defines the 

mated stellar mass ratio M^M,- 019(|Abt et al.[ 199 1|) Using outer northern boundary of the EON and separates M 43 from 

bispectrum speckle interferometry,| Preibisch et al. |(H999j) found the Huygens region . As inferred f r om me optical p olarization 

a third companion at ~ 470 mas with a flux ratio in the K-band map of the M 43 region obta i ne d by iKhallesse et alJ WM> . this 

of 0.03 + 0.02, implying a stellar mass ratio M 3 /M 1)2 ~ 0.07. dust lane is not a f oregroU nd obscuration but a wall of dense 

Although the less massive components may affect the spectrum material between M 42 and M 43. This wall of material ensures 

that it is quite unlikely that ionizing light from the Trapezium 

Send offprint requests to: ssimon@iac.es cluster stars reaches the nebular material ionized by HD 37061; 

* The INT and WHT telescopes are operated on the island of La 

Palma by the RGO in the Spanish Roque de los Muchachos Observatory 1 The Ha surface brightness in the EON is ~ two orders of magnitude 

of the Instituto de Astrofisica de Canarias. lower than the brightest regions in M 42. 
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Fig. 1. Color composite (RGB) image of the M42+M43 re- 
gion obtained from a combination of narrow band images taken 
with the WFC-INT (only CCDs #1, #2, and #4 are shown). 
The following color code was used: Ha (red), [S n] AA61 16+30 
(green), [Oiii]/15007 (blue). The field-of-view of the image is 
~ 22 arcmin x 33 arcmin. M43 is the roundish nebulosity at the 
north-east of the center. 



however, owing to the proximity of M 43 to the assumed north- 
ern boundary of the EON, a scattered light contribution (similar 
to the one described above) may still affect affect the nebular 
emission received from M43. 

Most studies of M43 have been performed at radio and 
far-infrared wavelengths. These studies were centered on the 
investigation of the global properties of t he nebula (i . e. ge- 
ometry, T R and n„. iMills & Shaved. Il968t iThumet all Il978t 
iThronson et al.L 119861: ISubrahmanvanl 1 19921) , and the charac- 
teristics of the dus t re-rad iation inside the nebula dThum et al.L 
Il978t ISmith et al.L Il987l) . As for optical images, radio con- 
tinuum observa t ions i ndicate that M43 has a circular shape 
(ISubrahmanvanLll992l). In addition, images at 60 /im pub lished 
by Thronson et al.l (1 19861) . and the kinematical study by lHanell 
allow one to associate the nebula with a shell-shaped cav- 
ity on the side of a dense molecular ridge. 



Its proximit^Q, relatively high surface brightness, simple ge- 
ometry, and isolated ionizing source make the Mairan's nebula 
an ideal object to investigate several topics of interest in the field 
of H ii regions and massive stars. To this aim, we collected a set 
of high quality observations comprising the optical spectrum of 
HD 37061, along with optical imaging and long-slit spatially 
resolved spectroscopy of M43. These observations were used 
to perform a comprehensive study of M 43 and its ionizing star, 
which is presented in two parts. In this first paper, we (1) obtain 
the stellar parameters of HD 37061, (2) investigate the presence 
of the scattered light in the region of M43 and its effect on 
the derived nebular properties, and (3) determine the physical 
conditions and nebular abundances in M 43 obtained at various 
distances to the center of the nebula. All this information is 
then used in a second paper (Simon-Diaz et al., Paper II, in 
prep.) to construct a specialized photoionization model of the 
nebula using as many observational constraints as possible. 
The main drivers of this second part of our study are to (1) 
test the reliability of the ionizing spectral energy distributions 
(SEDs) provided by the modern stellar atmosphere codes using 
nebular constra i nts, (2) investigate the effect of stellar pumping 
(iFerlandl [l999t iLuridiana et all 120091) on the nebular emission 
arising from the inner part of this H n region, and (3) test 
our understanding of the nebular energy budget temperature 
distribution. 

This paper is structured as follows. The observational 
data set is presented in Sect. [2] A quantitative spectroscopic anal- 
ysis of the optical spectrum of HD 37061 is performed in Sect. 
[3] The morphological characteristics of the nebula, along with 
its physical conditions and nebular abundances are determined 
in Sects. [4] and E]by analyzing of the Ha, H/3, [Oiii], and [Sn] 
narrow-band images, and the nebular spectra, respectively. Both 
nebular analyses also allow us to find several indications of an 
extended nebular emission, not directly related to M43, which 
must be adequately subtracted for a correct interpretation of the 
nebular images and spectroscopy. Finally, we outline the main 
results of this study in Sect. [6] 



2. The observational data set 

2.1. Nebular imaging 

The Wide Field CamercQ (WFC) attached to the Isaac Newton 
Telescope (INT) at the Roque de los Muchachos Observatory 
(La Palma, Spain) was used to obtain several narrow-band filter 
images of M 43 in 2004 February 1 6, and 2006 October 2 1 . The 
WFC-INT consists of four thinned EEV 2kx4k CCDs located 
at the primary focus of the telescope. Each CCD yields a field 
of view of 11.2' x 22.4', which results in a edge to edge limit 
of the mosaic (neglecting the ~1 arcmin inter-chip spacing) of 



2 The distance to HD 37061 and M43 can be assumed to be sim- 
ilar to that of its companion nebula M42. T wo different studies 
using Very Long Bas eline Array observations dMenten et all 120071 ; 
ISandstrom et al.LL2007l) determined the distance to some stellar objects 
in M42 to be 414±7pc and 389^} pc, respectively. These new mea- 
surements are ~ 10% lower than the 450 pc distance ofte n assumed for 
the Or ion nebular cloud (ONC). We refer the reader to iMenten et all 
d2007t) for a discussion of the diff erent dis t ances that have been deter- 
mined using different methods (viz. [Jeffries ( 2007 ): 392 pc, bas ed on the 
statistical properties of rotation in pre-MS- stars : lGoudi si dl982l) : 480 pc, 
for a long time considered the canonical distance). 

3 http : / /www . ing . iac . es/Astronomy/ instruments/wf c/ 
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M43 - Ha 




Fig. 2. H a, [S n] AA61 16+30, and [O m] A5QQ1 WFC-INT images of M 43. The size and position of the apertures extracted from the 
ISIS-WHT long-slit spectroscopic observations are also indicated. 



Table 1. Summary of the observations collected for the study. 



Nebular imaging 


Instrument 


Filter 


# lo (A) 


FWHM (A) 


(Telescope) 








WFC 


Ho- 


197 6568 


95 


(INT) 


Ho- redsh. 


228 6805 


93 


(0.33 arcsec/pix) 


H/? narrow 


225 4861 


30 




Hfi broad 


224 4861 


170 




[Om] 


196 5008 


100 




[S ii] 


212 6725 


80 


Nebular spectroscopy 




Grid 


A (A) Spect. range (A) 


Al (A/pix) 


ISIS 


R600B 


4368 3386-5102 


0.45 


(WHT) 


R600R 


6718 6064-7585 


0.49 


Stellar spectroscopy 




Grid 


lo (A) Spect. range (A) 


Al (A/pix) 


IDS 


H 2400B 


4320 4060 - 4590 


0.24 


(INT) 


H 2400B 


4800 4550 - 5070 


0.24 




H 1800V 


6400 6090 - 6760 


0.31 



34.2 arcmin. The pixel size is 13.5 /mi corresponding to 0.33 
arc sec/pixel. 

We centered M 43 on chip#4 of the WFC, and obtained im- 
ages of the region using the narrow-band filters H/3, [O in] /15007, 
Ha, and [S n] /L16716+30. The log of the observations, along 
with the characteristics of the various filters can be found in 
Table [T] While the [Om] and [Sn] images were only obtained 
to provide some qualitative information about the nebula, we 
wished to produce pure emission-line, flux-calibrated images for 
Ha and H/3. We hence obtained these images during a photomet- 
ric night, and included the Ha redshifted, and H/3 broad filters 
in our set of observations, to correct the Ha and H/3 images for 
the adjacent continuum. The exposure times considered for the 
images in the various filters were 4x60 s for Ha, 4x120 s for Ha 
continuum, 3 x 120s for H/5, 3x180 s for H/3 continuum, 5x12 s 
for [O m], and 5x12 s for [S n]. 

The reduction of the images was performed following 
the standard procedures (trimming, bias subtraction, flat- 
fielding, and alignment) with IRAfQ. Ha and H/? images 
were flux calibrated using observa tions of the spectrophoto- 
metric standard star BD+28°4211 (lLandolt & Uomotol 120071) 
at different airmasses. We obtained the continuum-subtracted 
Ha and H0 images following the procedure described in 
lLopez-Sanchez & Estebanl (120081) . 



4 IRAF is distributed by NOAO which is operated by AURA Inc., 
under cooperative agreement with NSF 
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Fig. 3. IDS-INT optical spectrum of HD 3706 1 . The H i and He i- 
ii lines used for the spectroscopic analysis of the star are indi- 
cated. 



The large field of view of the WFC also allowed us to include 
M42 in the images. A color-composite image of the M43+M42 
region obtained from a combination of the [Om] (blue), Ha 
(green), and [S n] (red) images is presented in Fig. [TJ Individual 
images in each one of these three filters, showing a square region 
of ~ 6 arcmin centered in M 43, are presented in Fig. |2] 



2.2. Stellar spectroscopy 

The spectroscopic observations of HD 37061 were carried out 
with the INT on 2006 August 30 and September 5. The 
Intermediate Dispersion Spectrograph (IDS) was used with the 
235 mm camera and two different gratings (see Table[TJi. We ob- 
served the spectral region between /1/14000 and 5050 A using the 
H2400B grating, which resulted in an effective resolving power 
R ~ 7500 (equivalent to a 0.24 A/pixel resolution and ~ 2.6 pixel 
fwhm arc lines). For the Ha region (~ /L16090-6760 A), the 
HI 800V grating was used, resulting in a similar spectral reso- 
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Table 2. Stellar parameters derived through the FASTWIND analysis 
of the optical spectrum of HD 37061. 



6500 7000 7500 

Wavelength (A) 

Fig. 4. ISIS-WHT optical spectrum of M43 (aperture #5). The 
Hi Balmer lines, along with Hei and other metal lines used in 
our study, are indicated. 



lution (0.31 A/pixel, R ~ 8000). With these configurations, three 
exposures were needed to cover the whole range. We obtained 
spectra with exposure times of 180s and 70s for the blue (2) 
and red (1) configurations, respectively, and obtained two spec- 
tra per spectral range. A large number of flat fields and arcs for 
the data reduction process were obtained. 

The reduction and normalization of the spectra was made 
following standard techniques, with IRAF and our own software 
developed in IDL. The signal-to-noise ratio of the reduced spec- 
tra is about 300-400, depending on the spectral range. Fig. [3] 
shows the global stellar spectrum, along with the H i and He i-n 
lines used for the quantitative spectroscopic analysis. 

2.3. Nebular spectroscopy 

A long-slit, intermediate-resolution spectrum of M43 was ob- 
tained on 2006 December 23 with the Intermediate dispersion 
Spectrograph and Imaging System (ISIS) spectrograph attached 
to the 4.2m William Herschel Telescope (WHT) at Roque de los 
Muchachos Observatory (La Palma, Spain). The 3.7' x 1.02" slit 
was located to the west of HD 37061, along the E-W direction 
(PA=90°, see Fig.©. Two different CCDs were used at the blue 
and red arms of the spectrograph: an EEV CCD with a configu- 
ration 4096 x 2048 pixels at 13.5 /an, and a RedPlus CCD with 
4096 x 2048 pixels at 15.0 yum, respectively. The dichroic used 
to separate the blue and red beams was centered on 5400 A. The 
gratings R600B and R600R were used for the blue and red ob- 
servations, respectively (see log of the observations in Table Q]). 
These gratings give a reciprocal dispersion of 33 A mm -1 in both 
cases, and effective spectral resolutions of 2.2 and 2.0 A, respec- 
tively. The blue spectra cover from /U3386 to 5102 A and the 
red ones from /L16064 to 7585 A. The spatial scales are 0.20" 
pixel -1 and 0.22" pixel -1 , respectively. The seeing during the 
observations was between ~0.5" and ~0.8". The exposure times 
were 3x300 s, in both the blue and red observations (obtained at 
the same time). 

The spectra were wavelength calibrated with a CuNe+CuAr 
lamp. The correction for atmospheric extinction was performed 
using the average curve for continuous atmospheric extinction 
at Roque de los Muchachos Observatory. The absolute flux 



r*(K) 


3 1000 ±500 


R(Ro) 


5.7: 


t0.8 


logg(dex) 


4.2 + 0.1 


logL/L G 


4.42: 


t0.12 


e(He) 


0.09 (assumed) 


M(M Q ) 


19 


±7 


logQ 


-15 (assumed) 


logQ(H°) 


47.2 


±0.2 



calibration was achieved by observations of the standard stars 
HD 19445, H600, and HD 84937. We used the standard IRAF 
TWODSPEC reduction package to perform bias correction, flat- 
fielding, cosmic -ray rejection, wavelength, and flux calibration. 
Fig.|4]shows an illustrative example of the global appearance of 
our nebular spectroscopic observations, also indicating the main 
nebular lines used in this study. 

3. Quantitative spectroscopic analysis of HD 37061 

The stellar parameters of the star were derived by visually com- 
paring the observed optical H i, He i, and He n line profiles with 
the synthetic ones resulting from t h e stell a r atmosphe r e cod e 
FASTWIND dSantolava-Rev et all Il997t iPuls et all l2005h . 
which is an established technique. Similar analys es, along with 
some notes on the methodolog y, can be found in iHerrero et al.l 
(I2002I) . iRepolust et all d2004l) . ISimon-Dfaz et al. I d2006l) . and 
references therein. 

To this aim, we constructed a grid of FASTWIND models 
with r eff , and logg, ranging from 28000 to 33000 K (500 K 
steps), and from 3.9 to 4.3 dex (0.1 dex steps), respectively. The 
microturbulence, the helium abundance, and the wind-strength 
parameteiQ were fixed to characteristic values for an early B 
dwarf star (f,=5kms -1 , e = 0.09, and log Q = -15, respectively) 
and metal abundanc es were assumed to b e solar (following the 
set of abundances bv lAsplund e7aTll2009l) . 

This technique requires the projected rotational velocity 
(v sin i) of the star to be previously established. All metal lines in 
the spectrum of HD 37061 are blended or very shallow because 
of the large v sin i o f the star; th erefore, we decided to apply the 
Fouri er method (c.f. lGravll976l see also lSimon-Dfaz & Herrerol 
|2007| for a recent application to OB-type stars) to the He i lines, 
and obtained a v sin i ~ 200 km s -1 . 

The synthetic H i and He i-n line profiles resulting from the 
models were then convolved with the corresponding instrumen- 
tal and rotational profiles and compared with the observations. 
The best fit was found for T e ff = 31000 K and log g- 4.2 dex 
(see Fig. |5j. Given the quality of the spectrum and the sensitivity 
of the used lines to the variation in r e ff and log g for this range 
of stellar parameters, an accuracy of 500 K and 0. 1 dex, respec- 
tively, could be achieved. As an example, in Fig.[5]we show the 
effect of a variation of +1000 K in the FASTWIND models. For 
this range of stellar parameters, the He n lines are decisive in 
constraining the T e $. 

We note the poor quality of the fits in the red wings of the H i 
and He i lines (but not in the Hen lines). These lines are likely 
affected by the spectrum of the cooler secondary spectroscopic 
component (see Sect.Q}. 

Once the stellar parameters have been determined, the stel- 
lar radius can be derived from the absolute visual magnitude 
(My) and the synthetic spectral energy distribution provided by 



Q = M/(Rv„y 
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Fig. 5. Fit of FASTWIND H i and He i-ii synthetic line profiles to 
the observed ones (red solid line, r e ff=31000 K, logg=4.2 dex). 
Blue dashed lines illustrate the effect of a variation of +1000 K 
in the FASTWIND models. Note that the red wings of the Hi 
and He i lines are affected by the cool companion. 



the stellar atmosphe re model. By using my =6. 84 and Ay-2.09 
dHillenbrandl Il997l) . and adopting a distance to the star of 
400+50 pc, we obtained an My = -3.3+0.3, quite consistent with 
the spectral type of HD 37061. The resulting stellar parameters, 
along with the total number of ionizing photons are summarized 
in Table [2] 



4. Qualitative and quantitative analysis of the 
nebular images 

Fig.[T]shows a color composite image of the M42+M43 region. 
The object under study, M43, appears as a roundish H n region 
centered around HD 37061, and well-separated from the Orion 
nebula (M42) by a dust lane (known as the northeast dark lane). 
The nebula has a diameter of ~4.5 arcmin (~ 0.65 pc at a distance 
of 400 pc), and is crossed by a dark lane oriented N-S in the 
eastern side (known as the M43 dark lane). This dark cloud is 
located in front of the nebula, blocking the nebular light coming 
from behind. 

Fig.|2]shows a closer view of the nebula, where three images 
of M 43, taken through filters isolating Ha (+[N n] /L16548+84), 
[S n] AA61 16+30, and [O m] /15007, are presented separately. 
The images show diffuse and extended emission beyond the lim- 
its of the nebula (mainly to the south and west). This emission 
can be more clearly seen in the [O ra] image. If the entire nebula 
were spherical with a density of 500 cm 3 , the size of the 2+ 
region would be ~ 25 % of the total size of the nebula^; however, 



6 Computed from a simple spherical Cloudy model using the ion- 
izing spectral energy distribution from a FASTWIND model with 
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Fig. 6. Three upper panels: Spatial distribution of Ha, 
[S n] AA61 16+30, and [Om] /15007 emission in a line passing by 
HD 37061 in the west-east direction (see Fig. O, obtained from 
the corresponding WFC-INT images. The dashed red line indi- 
cates the position of HD 37061. The vertical dotted line indicate 
the western limit of M43. Bottom panel: Ha/H/3 spatial distri- 
bution in the same line direction. Some dust features outlined by 
ISmith et alJ dl987h and the position and size of the ISIS-WHT 
slit are indicated. Horizontal dotted lines show the Ha/H/J value 
in the M43 west region. 



Fig. |2] shows that the [O ra] emission originates in a more ex- 
tended region (i.e. is unlikely to be relateqj to nebular material 
ionized by HD 37061). 

To more clearly illustrate this, we plot in Fig. |6](three upper 
panels) the spatial distribution of Ha, [Sn], and [Om] emission 
in a line passing by HD 37061 along the east-west direction, ob- 
tained from the corresponding WFC-INT images. Since we are 
only interested in the relative spatial changes, the three distri- 
butions are normalized to their maximum value. In all cases, a 
faint emission, not expected to be related to M43, is clearly de- 
tected beyond the limits of the nebula. In the [O ra] panel, we 
also superimpose the spatial distribution derived from the Ha- 
continuum image to show that the extended emission observed 
in the [O ra] image is not associated with the nebular continuum. 

4.1. Extinction 

The flux-calibrated, continuum-subtracted Ha and H/3 images 
(i.e. pure emission) were used to obtain the Ha/H/? flux ratio 
distribution along the diameter in the east-west direction passing 
by the ionizing star of M43 (see bottom panel in Fig.|6]l. Some 



l cff : 



:31000K, and log g = 4.2. 



7 This extended emission is more likely to be associat ed with the scat- 
tered n ebul ar light arising from the E ON, described by [O 'Dell & Gossl 
f2009h and lOTtel l & Harris (20l3). See a more detailed discussion 
about this possibility in Sect. l5.9l 
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interesting information about the extinction in the studied region 
can be extracted by inspecting of this figure: 

a) There is a jump in the Ha/H/3 flux ratio 140 arcsec west 
of HD37061 (i.e. western boundary of M43), indicating a 
clear difference in the type of emission arising from outside 
the limits of the nebula. We note that the lower Ha/H/3 ratio 
found in the "extended emission" region agrees with the hy- 
pothesis that this emission is associated with scattered light 
(see also Sect. 15.91 ), and is not necessarily an indication of a 
lower extinction in this region. 

b) The Ha/H/3 ratio in the western region of M 43 is fairly con- 
stant, though there is a small decrease close to the star. On 
the other hand, this ratio is larger in the eastern part. 

c) The M43 dark lane location (labeled as a foreground dust 
lane in Fig. [6j is clearly represented in the Ha/H/3 distribu- 
tion. 



ISmith et ail d 19871) presented a general view of the spatial 
distribution of dust in the M43 region. They found three regions 
where dust reradiation in 60 fim is concentrated (see a schematic 
map of dust features in their Fig l.c). The first one is associated 
with the M 43 dark lane (they called it the foreground dust lane), 
the second one (M43 east) is located to the east of HD 37061, 
between the ionizing star and the M43 dark lane, and a third 
one (M 43 front + M 43 west) follows the border of the nebula 
from the north to the west. This third feature of dust emission is 
separated ~ 60-70 arcsec from the star, leaving a "dust-empty" 
space in-between. We indicate in the bottom panel of Fig.|6]the 
location of these dust features using the same nomenclature as 
ISmith et al] 

The observed behavior of the Ha/H/3 flux ratio within M 43 
along the west-east directi on is perfectly corre lated with the dust 
distribution presented by ISmith et all (119871) . It is remarkable 
that, in contrast to expectations, the amount of extinction indi- 
cated by the Ha/H/3 flux ratio in the foreground dust lane is very 
low. The explanation is simple. This dark feature obscures the 
emission from M43, hence the observed light originates in the 
material located in front. This emission follows the trend indi- 
cated by the region marked as "extended emission". 

4.2. Total nebular Ha luminosity and surface brightness 

To derive the total nebular Ha luminosity we proceeded as fol- 
lows. First, the Ha image was corrected for the effects of dis- 
tance, extinction, and [N n] contamination^ Denoting with Fyi a 
the nebular flux received at Earth, the corresponding extinction- 
corrected Ha luminosity can be derived from: 



-Ha 



And 2 F»„lO c ^ ma \ 



(1) 



where d is the distance to the nebula, c(H/3) is the reddening 
coefficient, and /(Ha) is the value of the extinction function in 
Ha relative to H/3. 

The analysis of the Ha/H/3 image (Section [4. lb showed that 
the extinction is not constant across the nebula. The optimal 
strategy to follow would be to obtain c(H/3) for every pixel from 
the Ha/H/3 images, correct the Ha image, and then integrate the 
whole nebula to obtain the total L(Ha). However, when we at- 
tempted to follow this strategy we found that it introduces many 
sources of uncertainty. In addition, the nebular region hidden by 
the dark lane cannot be corrected in this way. We therefore de- 
cided to follow a different approach, taking into account what we 
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Fig. 7. Top and middle panels: integrated Ha luminosity profile 
and extinction corrected Ha surface brightness profile, respec- 
tively, obtained from the nebular image. Bottom panel: redden- 
ing coefficient obtained from the empirical analysis of the neb- 
ular spectra (see Sect. |5). In the upper panels we compare the 
L(Ha) and S(Ha) profiles of the original images (dotted lines) 
with those of images corrected for the contribution of the ex- 
tended emission (solid lines). Some data obtained from the ISIS- 
WHT spectroscopic observations are also included as triangles 
and diamonds (see text for explanation). 



learned from the inspection of the nebular images, and making 
use of the information extracted from the spectroscopic observa- 
tions (see Sect. [5)- We decided to integrate only the west quad- 
rant of the nebula, and multiply the resulting value by four (i.e. 
we assume the nebula is symmetric). Extinction in this quadrant 
is lower and more constant than in the other three. Furthermore, 
we have information about c(H/3) (~ 0.76) and the nitrogen con- 
tribution to the Ha image ([Nn] /i/l6548+84/Ha~ 0.59) for this 
region from the ISIS-WHT spectroscopy (see Sect. [5}- 

Once the distance, extinction, and nitrogen corrections had 
been applied to the continuum-subtracted, flux-calibrated Ha 
image, we obtained the L(Ha) inside circles of increasing sizes^ 
(centered on the star), and the surface brightness distribution 
within concentric rings of increasing radii. Fig. Qshows the mea- 
sured quantities as a function of distance to the star. Dotted lines 
show the initial determinations where the images were not cor- 
rected for the extended emission contribution. The Ha surface 
brightness distribution again shows that outside the limits of the 
nebula (9 > 160 arcsec) there remains a non-negligible, more or 
less constant emission. This implies that L(Ha) increases contin- 
uously, even outside the nebula (top panel in Fig. [7J. We there- 



Note that what we refer 
"Hff+[Nn] X16548+84" image. 



to as a "Ha image" is actually 



9 Actually, we integrated the mentioned quadrant and multiply the 
obtained value by four. 
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fore assumed a constant contribution of the extended emission 
to the surface brightness (~ 0.06 erg s -1 cm -2 , see middle panel 
of Fig. 13 and subtracted it from the Ha image. The resulting 
L(Ha) and S(Ha) distributions, calculated from the corrected 
images in the same way as described above, are indicated as solid 
lines in the top and middle panels of Fig. [7] These quantities 
behave as expected, i.e., the surface brightness becomes zero, 
and L(Ha) remains constant for R>R ne b- The resulting total Ha 
luminosity is (3.3+l.l)x 10 35 ergs -1 , where the associated un- 
certainty was calculated by taking into account uncertainties of 
12.5%, 10%, and 20% in distance, F(Har), and c(Hy8), respec- 
tively. 

The bottom panel of Fig.[7]shows the c(H/3) values obtained 
from the spectroscopic analysis of the 9 apertures extracted from 
the long-slit ISIS-WHT observations (see Sect. 15.2b . As for the 
Ha/Hfi images (bottom panel in Pig. [6j, we found that c(H/?) is 
fairly constant for 9 > 50 arcsec, but is somewhat smaller in the 
inner region (close to the star). We also include some informa- 
tion obtained from the spectroscopy in the middle panel of Fig. [7] 
(surface brightness profile). Red triangles correspond to the cal- 
culated S(Ha) in each aperture assuming a constant c(H/3) of 
0.76. The derived values follow the S(Ha) profile obtained from 
the Ha image. The middle panel in Fig. [7] also shows the corre- 
sponding S(Ha) values computed from the actual c(H/3) values 
derived from the spectroscopic analysis (black diamonds). These 
values deviate from the constant extinction-corrected values in 
the inner part of the nebula. We can conclude that the increase 
in S(Ha) in the inner part of the nebula is not necessarily caused 
by an increase in the emission measure or a deviation from the 
spherical geometry, but is more likely to be an effect of a non- 
constant extin ction or stellar pumping effects on_H recombina- 
tion lines (see iFerlandl 119991: iLuridiana etaii 120091) . The total 
L(Ha) we derived above must hence be corrected for this effect, 
resulting in (3.0+1.1) x 10 35 ergs -1 (about 10% below the pre- 
viously derived value). 

In Sect. [3] we derived the total number of ionizing pho- 
tons emitted by the central star (Q(H°) = 10 47 2±0 - 2 photons s -1 ). 
This would imply a maximum nebular Ha luminosity of 
(2.5+1.0)x 10 35 ergs -1 , a value that is in quite good agreement 
with the total Ha luminosity obtained from the analysis of the 
nebular images. This result is compatible to first order with M 43 
being a (mostly) ionization-bounded nebula in which dust does 
not absorb any significant fraction of the Lyman continuum pho- 
tons. 

We note that the stellar and nebular results would be incom- 
patible if we considered the total nebular Ha luminosity obtained 
from the images that are not corrected for the extended emission 
contribution. The total number of ionizing photons necessary to 
explain the derived nebular Ha luminosity would be at least dou- 
ble the value obtained from the analysis of the star. 

5. Empirical analysis of the nebular spectra 

5.1. More indications of extended nebular emission 

Fig. [8] shows the spatial distribution of the nebular emission 
in four selected lines (Ha, [O n] /U3720+30, He 1,16678, and 
[O m] /15007) and their respective adjacent nebular continua. 
These distributions were obtained from the ISIS@WHT spec- 
troscopic observations. The colored lines represent the measured 

10 We note that the Ho- surface brightness of the extended emission 
is not necessarily constant across the nebula; however, its variation is 
negligible compared to the Ho- emission from M43, hence can be con- 
sidered to be constant for the sake of simplicity. 
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Fig. 8. Spatial distribution of the nebular emission in four se- 
lected nebular lines (Ha, [O n] ^3720+30, HeiA6678, and 
[O m] /15007) and their respective adjacent nebular continua, 
obtained from the ISIS-WHT spectroscopic observations. The 
colored lines represent the measured line fluxes (blue), the ad- 
jacent continuum emission (red), and the continuum-subtracted 
line fluxes (black). The apertures used to extract nebular spectra 
at different distances from the central star are also indicated (see 
Sect. l5.21 >. Apertures labeled as vl, v2, v3, and v4, are those used 
to correct the other nine apertures for the contamination by the 
extended emission component. 



line fluxes (blue), the adjacent continuum emission (red), and the 
continuum-subtracted line fluxes (black). These figures display 
non-negligible emission beyond the limits of the nebula {9 > 150 
arcsec), that is not associated with the continuum emission. In 
most cases, this external emission remains fairly constant in the 
studied region (e.g. Ha, [O n] AA3720+30, Hei/16678, as well 
as [N n], [S n], and [S m] lines). However, in some cases (e.g. 
[O in]), this emission increases with increasing distance from the 
star. 

5.2. Aperture selection and line flux measurements 

We obtained one-dimensional (ID) spectra of regions of the neb- 
ula at different distances from the central star by dividing the 
long-slit used for the ISIS-WHT nebular observations into nine 
small apertures within the limits of the nebula (9 <150 arcsec) 
with the aim of obtaining: a) the radial distribution of the phys- 
ical properties of the nebula, b) information about nebular ionic 
and total abundances, and c) the radial distribution of critical 
nebular line ratios used as observational constraints on photoion- 
ization models of the nebula. 
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The size and position of the apertures are summarized in 
Tableland shown in Figs.|2]and[8] It is crucial to correct these 
spectra for the contribution of the extended emission compo- 
nent to obtain information about M43 itself. To achieve this 
aim, spectra from two extra-apertures outside the nebula were 
extracted to correct the line flux measurements of the apertures 
for the extended emission component. The size of the extra- 
apertures was selected to be the same as the corresponding aper- 
tures inside the nebula (11 and 16.5 arcsec for the three inner- 
most and the other outer apertures, respectively). Two different 
locations of these extra-apertures were considered depending on 
whether the considered nebular emission line had been emitted 
in the whole nebula (i.e. Hi, [O n], [N rx], [S n] and [S m] lines; 
apertures vl and v2 in Fig, [SJ, or only in the inner region (i.e. 
He i, [O m], and [Ar ni] lines; apertures v3 and v4 in Fig. [8}. For 
the latter set of lines, it is impossible to use the same apertures 
(vl and v2) as for the other lines, since the extended emission 
in these cases is not constant, but increases far away from the 
star. We hence considered two apertures as being close to the 
region where we were confident that the contribution of the neb- 
ular emission associated with M 43 had vanished. We expect the 
spectrum from these apertures to represent the properties of the 
extended emission close to the region where M 43 is emitting in 
the given nebular line. 

The extraction of ID spectra from each aperture was per- 
formed using the IRAF task apall. The same zone and spa- 
tial coverage was considered in the blue and red spectroscopic 
ranges. 

We detected H i and He i optical recombination lines, along 
with collisionally excited lines (CELs) of several ions, such as 
[O i], [O n], [O in], [N n], [S n], [S m], [Ne m], and [Ar m]. 
Line fluxes were measured using the SPLOT routine of the IRAF 
package by integrating all the flux included in the line profile be- 
tween two given limits and over a local continuum estimated by 
eye. Although we detected two [Ne m] lines, these were unsuit- 
able for analysis because they were extremely faint (in the case 
of /13869) or were severely blended with H 1/13967. In the case 
of [O n] /17330 A, we corrected the flux for the overlapping tel- 
luric emission line OH 8-3 P2 /17329.148 A by estimating its 
intensity from the observed OH 8-3 P2 T7329.148 A/OH 8-3 
PI /17340.885 A line ratio in standard star exposures. 

Each emission line in the spectra from the nine apertures was 
corrected for the extended emission component by subtracting 
the corresponding emission measured in the appropriate extra- 
aperture (vl to v4). Line intensities were then normalized to a 
particular H i recombination line present in each wavelength in- 
terval (H/3 and Ha for the blue and red spectra, respectively). To 
produce a final homogeneous set of line intensity ratios, the red 
spectra were then rescaled to H/3 using the theoretical Ha/H/3 
ratio I(Ha)/I(HjS) = 2.92 obtained for the physical conditions of 
T e = 7500 K and n e = 500 cnr 3 (see Sect. [53T >. 

5.3. Extinction correction 



We assume the extinction law derived by iBlagrave et al.l (120071) 
for the Orion nebula. Both Hn regions have a deficiency in 
small particles th at produces the relatively grey extinct i on ob - 
served in Or ion dMagazzu et all Il986t iBaldwin etaTT. Il99lh . 
Furthermore, iRodrfguezI (11999112002b did not find large extinc- 
tion variations in her optical spectroscopic studies of the M42 
and M 43 H ii regions. 

The reddening coefficient, c(H/3) was obtained by fitting the 
observed HS/H/3 and Hy/H/3 line intensity ratios - the three lines 



lie in th e same spectral range- to the theoretical ones com- 
puted bv lStorev & Hummed (1 19951) for T„ = 7500 K and n, — 500 
cm" 3 . As pointed out bv | Mesa-Delgado et al.l d2009albh . the use 
of a iBlagrave et all d2007l) law instead of the classical one by 
ICostero & Peimbertl (1 1970b produces slightly higher c(HyS) val- 
ues and also slightly different dereddened line fluxes depending 
on the spectral range. 

The final set of line intensities corrected for both the ex- 
tended emission component and the extinction are indicated in 
Tabled 

5.4. Uncertainties 

Several sources of uncertainties must be taken into account to 
obtain the errors associated with the line intensity ratios. We es- 
timated that the uncertaintjO in the line intensity measurement 
due to the signal-to-noise ratio of the spectra and the placement 
of the local continuum is typically 

~ 2% for F{A)/F(H/3)>0.5, 

~ 5% for 0.1< F(A)/F(H/3) <0.5, 

~ 10% for 0.05 < F(A)/F(H(3) <0.1, 

~ 20% for 0.01 < F(A)/F(Hf5)<0.05, 

~ 30% for 0.005< F(A)/F(H/3)<0M, and 

~ 40% for 0.001 < F(A)/F(H/3)<0.005. 

We did not consider any lines of weaker intensity than 0.001 x 
F(H/3). We note that then uncertainties indicated in Table prefer 
inly to this type of errors. 

By comparing the resulting flux-calibrated spectra of one 
of our standard stars with the corresponding tabulated flux, we 
could estimate that the line ratio uncertainties associated with 
the flux calibration is ~ 3% when the wavelengths are separated 
by 500- 1500 A and ~5% if they are separated by more than 
that. Where the corresponding lines are separated by less than 
500 A, the uncertainty in the line ratio due to uncertainties in the 
flux calibration is negligible. 

The uncertainty associated with extinction correction was 
computed by error propagation. The contribution of this uncer- 
tainty in the total error is again negligible when line ratios of 
close-by lines are considered (e.g. [S n] /16716/[S n] /16730). The 
final errors in the line intensity ratios used to derive the physical 
properties of the nebula were computed by quadratically adding 
these three sources of uncertainty. 

5.5. Physical conditions 

The electron temperature (T e ) and density (n e ) of the ionized 
gas were derived from classical CEL ratios, using the I RAF 
task temden of the package nebular (IShaw & Dufourl Il995h with 
updat ed atomic data for several ions (see iGarcfa-Roias et all 
l2009h . We computed n e from the [S n] A6717/i6731 and [O n] 
/13729//13726 line ratios and T e from the nebular to auroral [O n] 
/U(7320+30)//U(3726+29) line raticQ 

We could not directly determine r e ([0 m]) because the [O m] 
/14363 auroral line was not detectedbecause of the low ionization 



11 Although uncertainties in the line flux measurements depend on the 
line flux instead of on F(A)/F(H/3), the similarity between the H/3 flux 
measured in each aperture (see Table [3} makes this approach reason- 
able. 

12 In the case of M 43, T c ([0 n]) is free of any contamination to [O n] 
/17325 by recombination because of the low O 2+ abundance. 
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Table 3. Line intensities corrected from foreground emission and extinction (H/3 = 100), and results from the empirical analysis of 
the nebular spectra corresponding to the nine extracted apertures^ 1 '. 



Aperture 









Al 


A 2 


A3 


A4 


A 5 


A6 


A7 


A 8 


A9 


Center position (arcsec) 


15.80 


26.80 


37.80 


51.55 


68.05 


84.55 


101.05 


117.55 


134.05 




Size (arcsec) 


11.0 


11.0 


11.0 


16.5 


16.5 


16.5 


16.5 


16.5 


16.5 


A (A) 


Ion 


Mult. 








7(i)//(H/5) 












3726.03 


[O ii] 


IF 


93.5+1.8 


99.7±2.0 


106±2 


104±2 


104±2 


108±2 


117±2 


129+3 


150+3 


3728.82 


[O ii] 


IF 


84.2±1.7 


91.5±1.8 


98.1±2.0 


99.5±2.0 


98.9±2.0 


103±2 


114+2 


120+2 


143+3 


3835.39 


Hi 


H9<2) 


7.50±0.75 


7.94±0.79 


8.28±0.83 


7.91±0.79 


7.83±0.78 


7.57±0.76 


7.98+0.80 


7.37+0.74 


7.42+0.74 


4068.60 


[S ii] 


IF 


1.58±0.32 


1.80±0.36 


1.86±0.37 


2.43±0.49 


2.70±0.54 


3.30±0.66 


3.84±0.77 


4.77+0.95 


5.63+0.56 


4076.35 


[S ii] 


IF 


0.64±0.19 


0.67±0.20 


0.76±0.23 


0.74±0.22 


0.71±0.21 


1.00±0.30 


1.16±0.23 


1.29+0.26 


1.24+0.25 


4101.74 


Hi 


HS 


26.2±1.3 


26.2±1.3 


26.0±1.3 


25.7±1.3 


25.6±1.3 


25.6+1.3 


25.7+1.3 


25.6+1.3 


25.8+1.3 


4340.47 


H l 




46.2±2.3 


46.2±2.3 


46.5±2.3 


46.9±2.3 


46.9±2.4 


47.0±2.4 


46.8±2.3 


46.9+2.4 


46.8+2.3 


4471.09 


He i 


14 


1.09±0.22 


1.35±0.27 


1.02±0.20 


0.19±0.08 


— 


— 


— 


— 


— 


4861.33 


Hi 


H/i 


100±2 


100±2 


100±2 


100±2 


100±2 


100+2 


100+2 


100+2 


100+2 


4958.91 


[O m] 


IF 


2.52±0.50 


0.98±0.29 


0.15±0.06 








0.93+0.28 


0.88+0.26 


7.59+0.76 


5006.94 


[O m] 


IF 


7.86±0.79 


2.78±0.56 


0.25±0.10 




0.35±0.14 


0.13±0.05 


2.78±0.56 


3.06+0.61 


25.42+1.27 


6300.30 


[O i] 


IF 


0.64±0.03 


0.64±0.03 


0.60±0.03 


0.86±0.04 


l.O5±0.05 


1.30±0.06 


2.47±0.12 


3.50+0.17 


5.39+0.27 


6312.10 


[S hi] 


3F 


0.80±0.04 


0.70±0.03 


0.69±0.03 


O.58±0.03 


0.47±0.02 


0.42±0.02 


0.45±0.02 


0.38+0.02 


0.36+0.02 


6548.03 


[N ii] 


IF 


36.31±1.82 


38.4±1.9 


41.3±2.1 


41.8±2.1 


42.0±2.1 


44.2±2.2 


45.7±2.3 


50.2+1.0 


54.1 + 1.1 


6562.82 


Hi 


Ha 


ZVZ. UU±J. 54 


zv_±d 


IQl+f. 
ZVZ+D 


IQl+f. 
ZVZ+D 


IQl+f. 


IQl+f. 
ZVZ+D 


IQl+f. 
ZVZ+D 


IQl+f. 
ZVZ+D 


IQl+f. 
ZVZ+D 


6583.41 


[N ii] 


IF 


1U/. 1Z+Z. 14 


1 1 1+1 
1 1 J+Z 


1 11+1 
1ZZ±Z 


1 11+1 
1ZJ+Z 


1 1A+1 
1Z4±Z 


i it\+i 


1 1^ + 1 


1 A Q + 1 


1 JV+J 


6678.15 


He i 


46 


1 . 1 1 ±u.uo 


1 . 1 J+U.UO 


U, BV+U.U^t 


n ic\+c\ ni 

U.ZU+U.U 1 








U.Uj+U.U 1 


U.ZO+U.U 1 


6716.47 


[S ii] 


2F 


18.79±0.94 


20.2±1.0 


22.2±1.1 


25.8±1.3 


29.4±1.5 


36.9+1.8 


42.3+2.1 


49.1+2.5 


58.8+1.2 


6730.85 


[S ii] 


2F 


20.02±1.00 


21.4+1.1 


23.0±1.2 


26.6±1.3 


30.4±1.5 


38.7+1.9 


42.8+2.1 


53.3+1.1 


60.8+1.2 


7065.28 


He i 


10 


O.92±0.05 


0.85±0.04 


O.65±0.03 


0.15±0.01 








0.06+0.01 


0.09+0.01 


0.34+0.02 


7135.78 


[Ar m 


] IF 


2.90±0.14 


2.66±0.13 


1.70±0.08 


0.31 ±0.02 


0.03±0.01 


0.03±0.01 


0.19+0.01 


0.24+0.01 


1.16+0.06 


7319.19 


[On] 


2F 


1.74±0.09 


1.67±0.08 


1.67±0.08 


1.51±0.08 


1.52±0.08 


1.66±0.08 


1.80+0.09 


2.53+0.13 


2.67+0.13 


7330 20 


[On] 


2F 


1.46±0.07 


1.41±0.07 


1.34±0.07 


1.24±0.06 


1.26±0.06 


1.39±0.07 


1.50+0.08 


2.07+0.10 


2.20+0.11 






com 


0.60±0.09 


0.49±0.08 


0.51±0.05 


0.75±0.06 


0.76±0.06 


0.79±0.07 


0.77+0.06 


0.73+0.06 


0.99+0.07 








1.508xl0~ 12 


1.443xl0~ 12 


1.185xl0~ 12 


1.861xlO~ 12 


1.770xl0~ 12 


1.765xl0~ 12 


1.271xl0~ 12 


1.150xl0~ 12 


7.131xl0~ 13 






nMO n]) 


560±50 


520±5O 


500±40 


440±40 


450±40 


450±40 


420+40 


510+50 


460+40 






"MS U]) 


650±170 


630±170 


570±160 


560±150 


560±150 


600±160 


520+150 


690+140 


580+60 






T c ([0 ii]) 


7850±160 


7600±150 


7360±140 


7260±130 


7270±140 


7420±140 


7440+140 


8070+180 


7810+180 




F e ([0 m]f> 


7360±1700 


7000±1650 


6660±1600 


















He + /H + 


10.44±0.02 


10.45±0.03 


10.34±0.03 


















0+/H+ 


8.34±0.05 


8.45±0.05 


8.55±0.05 


8.58±0.05 


8.58±0.05 


8.54±0.05 


8.57+0.05 


8.42+0.05 


8.56+0.06 






2+ /H + 


6.97±0.20 


6.64±0.22 


5.83±0.26 


















N + /H + 


7.63±0.04 


7.70±0.04 


7.78±0.04 


7.80±0.04 


7.80±0.04 


7.80±0.04 


7.81+0.04 


7.73+0.04 


7.81+0.04 






S + /H + 


6.28±0.04 


6.35±0.04 


6.43±0.04 


6.51±0.04 


6.56±0.04 


6.64±0.04 


6.68+0.04 


6.66+0.04 


6.77+0.04 






S 2+ /H + 


6.78±0.04 


6.80±0.04 


6.87±0.04 


6.83±0.04 


6.73±0.04 


6.64±0.04 


6.66+0.04 


6.39+0.05 


6.44+0.05 






Ar 2+ /H + 


5.79±0.16 


5.82±0.17 


5.69±0.17 


















O/H 


8.36±0.05 


8.45±0.05 


8.55±0.05 


8.58±0.05 


8.58±0.05 


8.54±0.05 


8.57+0.05 


8.42+0.05 


8.56+0.06 






N/H 


7.65±0.08 


7.71 ±0.08 


7.78±0.08 


7.80±0.04 


7.80±0.04 


7.80±0.04 


7.81+0.04 


7.73+0.04 


7.81+0.04 






S/H 


6.90±0.04 


6.93±0.03 


7.01±0.03 


7.00±0.03 


6.96±0.03 


6.94±0.03 


6.97+0.03 


6.85+0.03 


6.94+0.03 



(1) The errors in the line fluxes refer only to uncertainties in the line measurements (see text). F(HJ3) in erg cm 2 s 1 ; w e in cm 3 ; T c in K; ionic abundances in log (X + '/H + ) +12. 

(2) Blended with He i ,13833.57 line. 

(3) r c ([0 hi]) obtained using the empirical relation between 7" c ([0 n]) and T c ([0 in]) obtained from the data by Garcia-Rojas & Esteban 1 2007) (see text). 



degree of the nebula. As an alternative, we used a empirical rela- 
tion between T e ([0 in]) a nd r R ([Q nl) obtained from the ne bular 
information published by iGarcfa-Roias & Esteban I (120071) for a 
sample of Galactic H n regions! 13 ! 

r e ([ODI]) = [r e ([OII])-(2640+1270)]/(0.7105 + 0.1524).(2) 

This f it is almost identica l to that obtained by iPilvugin et al.l 
j2006l) and lGarnettl i 19921) for integrated spectra of giant extra- 
galactic H ii regions and photoionization models of giant extra- 
galactic Hn regions, respectively. We only computed r e ([0 m]) 
for the three innermost apertures. For more distant apertures, 
the 2+ emission (if any) is not associated to M43 itself, but 
to the residual emission of the diffuse component (see Sects. [4] 
andliUTi. 

The methodology used to determine the physical conditions 
was as follows: we assumed a representative initial value of T e 

13 The data for NGC3603, a giant Galactic Hn region, was discarded 
to obtain this relation. 



of 10000 K and computed the electron densities. The value of 
n e was then used to compute r e ([0 n]) from the observed line 
ratios, and T e ([0 m]) from equation 2. We iterated until conver- 
gence to obtain the final values of « e and T e . Uncertainties were 
computed by error propagation. The final « e ([0 n]), n e ([S n]), 
r e ([0 ii]), T e ([0 in]) estimate, along with their uncertainties are 
indicated in Table [3] 

In general, densities derived from the [O n] line ratio are 
about 100 cirT 3 lower than those derived from [S n] lines, but 
consistent within the errors. 

5.6. Chemical abundances 

To derive He + /H + , we used two observed lines of He i at /L1447 1 
and 6678. Case B emissivities were taken from the colli sion 
less (low-density limit) calculations bv lBauman et al. I J2001 us- 
ing an on line available codfQ The collisional-to-recombination 
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contribution was estimated from lKingdon & Ferlandl ( 1995|). us - 



ing the interpolation formula provided by Por ter et al. ( 20071) . 
The effective recombinati on coefficients for H + were taken from 
IStorev& Hummed (I 19951) 

Ionic abundances of N + , + , 2+ , S + , S 2+ and Ar 2+ were 
derived from CELs, using the IRAF task ionic of the package 
nebular. We assumed a two-zone scheme, in which we adopted 
r e ([0 n]) for the low ionization potential ions (N + , + , S + and 
S 2+ ) and the value of r e ([0 m]) derived from Eq.[2]for the high 
ionization potential ions (0 2+ and Ar 2+ ). In both cases, we as- 
sumed n e ([0 n]) instead of n e ([S n]) because the first indicator 
is more representative of the whole nebula, while n e {[S n]) pro- 
vides insight mostly into the density of the material close to the 
ionization edge of the nebula. For 2+ , Ar 2+ , and He + , we only 
have available lines in the three innermost apertures. Owing to 
the low T e s of the ionizing star, these relatively high ionization 
species are only present in regions very close to the star. 

We then derived the total abundances of O, N, and S for each 
aperturePl As a consequence of the low ionization degree of the 
nebula, M 43 has a very small amount of 2+ in comparison with 
the dominant species, + , except in the innermost zones where a 
non-negligible fraction of 2+ is present. We note that for the in- 
termediate apertures (apertures 4 to 6), the O abundance is given 
by + /H + . In the case of S, the abundance of this element is 
given by the sum of S + and S 2+ , because S +3 is not expected to 
be present in this nebula. Similarly to 2+ , N 2+ is only expected 
to be present in the innermost parts of the nebula. For the cor- 
responding apertures (1-3), we applied the classical ionization 
correction factor (ICF) scheme of N (i.e. N + /0 + =N/0) to correct 
for the N 2+ abundance (note that this correction is very small). 
As for O, for the intermediate apertures the N abundance is di- 
rectly obtained from the N + abundance. 

The ionic abundances of He + + , 2+ , N + , S + , S 2+ , and 
Ar 2+ , along with the total abundances of O, N, and S are shown 
in Table [3] We stress again that in this nebula the total abun- 
dances of O, N, and S can be obtained directly from observable 
ions (no ICFs are needed) from the intermediate apertures. 



5. 7. Comparison of O, N, and S abundances derived in M 43 
and the Orion nebula 

Given that the ionized gas from both nebulae, M42 and 
M43, comes from the same molecular cloud, it is reason- 
able to assume that both have very similar total abundances. 
In this section, we compare the O, S, and N abundances 
obtained from our study of M 43 with those determined by 
iGarcfa-Roias & Estebanl (120071 GRE07. an upd ate of Esteban 
et al. 2004) and ISimon-Dfaz & Stasihskal d2010l SDS10) from 
the analysis of the nebular spectrum of the Orion nebula. Table|4] 
summarizes the final abundances indicated in the three studies. 
We note that GRE07, and SDS10 present analyses of the same 
spectroscopic dataset, but use different atomic data and ICFs. 
Our study of M 43 uses the same atomic data as that considered 
by GRE07. 



5.7.1. Oxygen 

For this element, where the total abundance is obtained directly 
from observable ions (O + and O + +0 2+ in the case of M 43 and 
M42, respectively), good agreement is found within the errors. 



Table 4. Summary of O, N, and S abundances (from CEL) resulting 
from the spectroscopic analysis of M 43 (this work) and the Orion neb- 
ula (GRE07, SDS10). For M43, we assumed the mean value of abun- 
dances derived from apertures 4-6 (the uncertainty given by the values 
indicated in Table[3]for each aperture). 



Element 


M43 


M42 (GRE07) 


M42 (SDS10) 


O 


8.57±0.05 


8.54+0.03 


8.52+0.01 


S 


6.97±0.03 


7.04+0.04 


6.87+0.04 


N 


7.80+0.04 


7.73+0.09 


7.90+0.09 



5.7.2. Nitrogen 

For this element, GRE07 and SDS10 evaluate different values. 
This discrepancy is caused by the assumed ICF(N + ). While 
GRE07 obtained this ICF from the widely used empirical re- 
lation N + /0 + =N/0, SDS10 determined the ICF from a model 
fitted to the data of M42 (including a detailed description of 
the ionizing SED of Q'OriC predicted by the stellar atmo- 
sphere code WM-basic, iPauldrach et~aT1,l2001l) . The GRE07 and 
SDS10 abundances differ by 0.17 dex. Interestingly, our deter- 
mination for M43, in which no ICF is needed produces an in- 
termediate value, and is in perfect agreement with the B -type 
stars abundance derived by iNieva & Simon-Dfazl (1201 ll) . This 
introduces a new impo rtant constraint on the nebular C EL/RL 
abundance conundrum dSimon-Dfaz & S tasihskal [20101) . 



5.7.3. Sulphur 

GRE07 and SDS10 derive different values (0.17 dex differ- 
ence) of the total S abundance in M42. In this case, the 
cause of the difference is not, however, the assumed ICF, 
but the use of different sets of atomic data. Both compu- 
tations considered the same data f or the collision strength 
dKeenan et all Il996t iTaval & Guptall 19991) . but different tran- 
sition probabi lities: while GRE07 based t heir calculation 



in data from iMendoza & Zeipp"en (Il982albl) fo r [S nl and 
[Sml, respectively, SDS10 use d iFroese Fischer] (|2004 and 



15 Owing to the low ionization degree of M43, the ionization correc- 
tion factors for He + and Ar 2+ are large and very uncertain. We did not 
compute the total abundances for these two elements. 



iFroese Fischer & Tachievi (120041) values. For M 43, we obtained 
12 + log(S/H)=6.97+0.04 using the same set of atomic data as 
GRE07. This value is slightly lower, but in agreement within the 
errors, than the determination by GRE07, and somewhat higher 
than the value provided by SDS10. The influence of the atomic 
data considered for this element is hence relevant. 

5.8. Comparison with previous studies of M43 

The only nebular abundance studie s foun d in the literature for 
M 43 are those o f iRodrfguez I d 1 9991 120021) . who analyzed long- 
slit optical spectroscopic data for five different slit positions 
along the nebula, deriving the physical conditions (T e and n e ), 
and O, S, CI, N, Ar, He, C, and Fe abundances of the nebular gas. 
lO'Dell & Harris! (1201 Oh analyzed a set of optical spectrophoto- 
metric observations of M 43 as part of a more general study of 
the EON, and derived the physical conditions of the nebula. We 
note that none of these observations were corrected for the dif- 
fuse emission component as in our study. 

For consistency, we reanalyzed the data of lRodrfguez I (1 19991) 
using the same atomic data as in our work. We obtained 
T e ([N ii]) between 7800 K and 80 00 K for slit positions simi- 
lar to ours. lO'Dell & Harrisl(l2010l) obtained very similar results 
(r e ([N n]) between 7780 K and 7940 K) from several aperture 
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Fig. 9. Upper panel: Difference between the non-extinction- 
corrected Hy/H/? and H5/H/3 ratios from the C and NC spectra 
(see text for explanation). Lower panel: He 1 7065/6678 line ratio 
from the NC spectra (9 apertures). The mean and standard devi- 
ation values corres ponding to the bright H uygens region (ob- 
tained from data bv lO'Dell & Harris! 120101) are also indicated as 
a dashed and dotted lines, respectively. 



extractions of a slit that crosses the nebula in the east-west di- 
rection. Both authors derived « e ([S n]) ~500-600 crrT 3 , which 
is very similar to the value derived in this work. 

We compared our abundances with thos e resulting f r om th e 
reanalysis[]3 of the observational dataset bv lRodrfguezl (119991) . 
The derived O abundances are 0.14-0.23 dex lower than those 
resulting from the analysis of our spectra. As we indicate in Sect. 
15.91 this effect is produced by the contamination of the M43 
emission by the diffuse component. 

Finally, we compared our estimated r e ([0 n]) with pre- 
vious determinations from radio continuum observations. 
ISubrahmanvanld 19921) derived r e =9000 + 1700 K by combining 
330 MHz and 10.7 GHz co ntinuum observations and a model 
of an isothermal H n region; iMills & Shaver! (fl%8l) estimated a 
value of r e ~8000 K from 408 MHz continuum observations; fi- 
nally, [Thurnetal] dl978l) derived a much lower r e ~6700 K from 
the HI 91a' radio recombination line to continuum ratio. In gen- 
eral, the different estimates, except for Thum et al., agree within 
the uncertainties. 



5.9. Scattered nebular light from the EON affecting the 
spectrum of M 437 

lO'Dell & Gosi (l2009t) and lO'Dell & Harrisl (1201 Oh demon- 
strated that scattered light from the bright Huygens Region of 
M42 can affect t he physical conditions d erived from emission 
fines in the EON. lO'Dell & Harris! (1201 Oh performed a detailed 
spectroscopic study of several zones of the EON, including the 
M 43 region. These authors could not find conclusive evidence 
of signatures of scattered light from the Huygens region affect- 
ing the nebular emission from M 43 (except for the case of an 
aperture located in the M43 dark lane). Our study (based on 
more detailed observational dataset) allowed us to observation- 
ally prove that there is an extended diffuse emission not pro- 



16 The observations by [O 'Dell & Harrisl (l2010h do not include 
[O n]/Li3727, 7325 lines, hence the determination of the oxygen abun- 
dance is not possible. 



duced by ionization of the surrounding material by NU Ori af- 
fecting the whole nebula. 

To investigate whether this extended emission relates to scat- 
tered light from the Huygens region, we compare in Fig.|9](upper 
panel) the non-extinction-corrected Hy/H/3 and H£/H/3 line ratios 
with (C) and without (NC) correction for diffuse emission for 
apertures from Al to A9. The H6/H/3 line ratio is more affected 
than Hy/H/?. This is consistent with the diffuse extended emis- 
sion being produced by dust scattered light, a well-know prop- 
erty of this type of emission is that blue lines are more affected 
than red lines. An additional argument supporting this hypothe- 
sis is presented in the bottom panel of Fig. [9] The value of the 
He 1 7065/6678 ratio, obtained from the NC spectra of M 43 out- 
side the He + Stromgen sphere associated with NU Ori agrees 
with the associ ated average value in th e Huygens region taken 
from the data of lO'Dell & Harrisl (fcOlOlP l 

We also investigated the effect that scattered light contami- 
nating the spectrum of M 43 has on line ratios used to determine 
the physical conditions and chemical abundances (as well as 
constraints of the photoionization models presented in Paper II). 
Fig. [10] shows several line ratios obtained from the C and NC 
spectra. Although the effect is negligible for the density con- 
straints and the S + and O + abundance indicators, the influence 
of scattered light on the ionization structure constraints and the 
He + , N + , Ar 2 +, S 2 +, and O 2+ abundance indicators is not neg- 
ligible (especially in th e outer parts of the nebul a). In addition, as 
already pointed out by|Q 'Dell & Harrisl (120101) . our comparison 
identifies an important effect on r e ([On]). If the scattered light 
contribution is not removed, the spectroscopic analysis overesti- 
mates the T e and hence underestimates of the abundances. 

To quantitatively illustrate this effect on the computed phys- 
ical conditions and abundances, we present in Table[5]the results 
of two empirical analyses of the C and NC spectra from aper- 
tures 4 and 5. We also compare the result of the analysis of the 
NC integrated spectrum (extracting the global long slit spectra). 
As pointed out before, n e is basically similar but the NC spec- 
tra result in a higher T e and, consequently lower abundances. In 
particular, ionic abundances can be affected by ~ 0.05 -0. 15 dex. 

An adequate correction of the scattered light contribution in 
the spectra of M 43 is hence important for a reliable empirical 
analysis of the nebular spectra of this H n region. 



6. Summary and conclusions 

M 43 is a close-by, bright H n region of simple geometry, ionized 
by a single star. These characteristics makes M43 an ideal ob- 
ject for investigating several topics of interest in the field of H n 
regions and massive stars. 

In a series of two papers, we present a combined, com- 
prehensive study of the nebula and its ionizing star using as 
many observational constraints as possible. In this first part 
of the study, we have introduced the observational dataset, 
obtained the stellar parameters of HD 37061, obtained useful 
information from the nebular images, and analyzed the nebular 
spectra extracted from apertures located at various distances 
from the central star. All this information is then used in Paper 
II to construct a customized photoionization model of the nebula. 

In this paper, we have found observational evidence of a 
diffuse and extended emission in the region where M43 is lo- 
cated that is not associated to material ionized by HD 37061. 



17 We considered extractions corresponding to slits 2 to 15 that are 
labeled as "inner". 
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Fig. 10. Observed line ratios obtained from the corrected (black diamonds) and non corrected (red squares) spectra from extended 
emission component. 



Table 5. Summary of physical conditions, ionic and total abundances 
derived for apertures 4 and 5 with and without correcting for extended 
emission component and for the longslit extraction in M43. (1> 







Uncorrected 


Corrected 


Uncorrected 






A4 


A5 


A4 


A5 


Total Slit 


« c (cm 3 ) 


[On] 


475 


485 


440 


450 


510 




[Sn] 


560 


565 


560 


560 


600 


r e (K) 


[On] 


7650 


7700 


7260 


7270 


7900 




[O in] 










7450 


6(X + ') 


o + 


8.45 


8.44 


8.58 


8.58 


8.38 




o ++ 










7.45 




N + 


7.71 


7.70 


7.80 


7.80 


7.66 




S + 


6.42 


6.46 


6.51 


6.56 


6.42 




S ++ 


6.77 


6.70 


6.83 


6.73 


6.69 




Ar ++ 










5.64 


e(X) 


O 


8.50 


8.49 


8.58 


8.58 


8.43 




N 


7.75 


7.75 


7.80 


7.80 


7.71 




S 


6.93 


6.90 


7.00 


6.96 


6.88 



7" e , ionic, and total abundance values are based on the assumption « c =« e ([0 n]) as 
discussed in Sect. 15.61 



Our nebular observations have allowe d us to ascertain the 
most likely origin of t his emission (e.g. lO'Dell & GossL 120091: 
lO'Dell & HarrisLfeOlOl and references therein), namely that light 
emitted in the Huygens region (the central, brightest part of the 
Orion nebula) is scattered by dust. We have also shown the im- 
portance of an adequate correction of this scattered light from 
the imagery and spectroscopic observations of M 43 to a proper 
determination of the total nebular Ha luminosity, the nebular 
physical conditions, and chemical abundances. In particular, we 
illustrated that an overestimate of T e by ~ 400-500 K, hence an 
underestimate of abundances by 0.05-0.10 dex, result from the 
empirical analysis of our spectroscopic dataset when the spectra 
are not corrected for the scattered light contribution. In addition, 
the derived total nebular Ha luminosity may be overestimated 
by a factor ~ 1.5. 

The quantitative analysis of the optical spectrum of 
HD 37061 with the stellar atmosphere code FASTWIND lead to 



r eff = 3 1000±500 K, log g = 4.2+0. 1 , and log Q(H°) = 47.2 + 0.2 
(assuming a distance of 400 pc). 

The analysis of the Ha and H/3 images indicate a non- 
constant extinction distribution within the nebula that is wel l 
correlated with the dust features indicated bv lSmith et al.l(ll987l) . 
Once the Ha image is corrected for extinction and diffuse emis- 
sion, a total nebular Ha luminosity of (3.0+1.1) x 10 35 ergs~' is 
obtained. This value is compatible with the ionizing flux of the 
star, implying that the nebula is (mostly) ionization bounded. 

We extracted nine apertures from a long-slit located to the 
west of HD 37061 (east-west direction) to obtain the spatial dis- 
tribution of the nebular physical conditions (temperature and 
density) and ionic abundances (He + , + , 2+ , N + , S + , S 2+ , 
Ar 2+ ). Since it is important to correct these spectra for the con- 
tribution of scattered light from the Huygens region, we also 
consider two apertures outside the nebula that are used for this 
aim. 

For sulfur, oxygen, and nitrogen, we have been able to de- 
termied total abundances directly from observable ions (no ICFs 
are needed). The derived abundances, 8.57 + 0.05, 6.97 + 0.03, 
and 7.80 + 0.04, respectively, are compared with previous deter- 
minations in the Orion nebula. Although an overall agreement 
is found, our study illustrates the importance of the atomic data 
and, specially in the case of M 42, the considered ICFs. 
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